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We investigate the effects of the next-nearest-neighbor (t') and the third- nearest-neighbor (t") 
hopping terms on superconductivity correlation in the 2D hole-doped extended t — J model based 
on the variational Monte-Carlo, mean-field calculation, and exact diagonalization method. Despite 
of the diversity of the methods employed, the results all point to a consistent conclusion: While 
the d— wave SC correlation is slightly suppressed by t' and t" in underdoped regions, it is greatly 
enhanced in the optimal and overdoped regions. The optimal T c is a result upon balance of these 
two opposite trends. 

PACS numbers: 74.20.-z 



Right after the discovery of high temperature super- 
conductors, the two-dimensional (2D) t — J model has 
been proposed to provide the mechanism of supercon- 
ductivity (SC) |l| . This idea quickly gained momentum 
when variational calculations showed that the doping de- 
pendence of pairing correlation |2j y| and the phase dia- 
gram of antiferromagnetic (AFM) phase and SC 0] seem 
to agree with experimental results fairly well. However, 
recently calculations beyond the variational method 
have challenged the notion that pure 2D t — J model 
without including other interactions is enough to explain 
the high values of T c . Although this issue is yet to 
be settled B, B, la, there are results by band-structure 
calculations |9l flQ| and experimental analysis ^lj that 
hopping beyond nearest neighbors is essential to raise 
T c . In fact they found highest T c<max for different mono- 
layer cuprates strongly correlates with t'/t, where t' is 
the second nearest-neighbor ho pping amplitude. This 
contradicts with previous results|l2ul3| of exact calcula- 
tions that for the hole doped systems, introducing t' into 
the t-J model will actually reduce pairing. Although the 
last results are studied for systems doped only with a few 
number of holes, all these conflicting results raised a very 
serious challenge to the t — J type models. 

There are many experimental and theoretical results 
to support the presence of t' and possibly also t" , the 
third nearest neighbor hopping, in cuprates. The topol- 
ogy of large Fermi surface seen by ARPES0, and 
the change of sign of Hall coefficient as a function of 
doping 16] can best be understood in the presence of t' 
and t". The single hole dispersion observed by ARPES 
and the difference between hole- and electron-doped 
systems also support the presence of these terms. 

In this letter, we will use variational approach sup- 
plemented by slave-boson mean-field (MF) calculations 
and exact diagonalization (ED) method to show that the 
presence of t' is indeed important for enhancing pairing 



beyond the underdoped regime. The largest values of 
pairing correlation obtained are proportional to t'/t up 
to t'/t = —0.3 ~ —0.4. In addition we will show that the 
decrease of pairing correlation at very large hole density 
is related to the change of Fermi surface topology. The 
conflicting results between theories and experiments dis- 
cussed above are naturally resolved within the extended 
t — J model. 

The Hamiltonian of the extended t — J model is 

H =H t + H. r = (1) 

ij <i,j> 

where ty = t, t' , and t" for sites i and j are nearest, next 
nearest, and the third nearest neighbors, respectively, 
and tij = for longer distance. Ci jCr = (1 — ni-a)c%,ai 
satisfies the no-double-occupancy constraint. In our no- 
tation for hole doped materials, t'/t is negative while 
t" = —t'/2 most of the time. 

The trial wave function (TWF) used in this study is 
the d x 2_ y 2 resonating- valence-bond wave function 

\*)=P G ]l(u k +v k cl A cl kil ) |0) (2) 
k 

with u k /v k = A fc /(e fe + e\ + A 2 .), e k = -2t(cosk x + 
cosky) — At' v cosk x cosk y — 2t"{cos2k x + cos2k y ) — /i, and 
Afc = 2A(cosk x — cosky). Here the projection opera- 
tor Pq enforces the constraint of one electron per site. 
In addition to A and fi, we have included two impor- 
tant variational parameters t' v and t v which determine 
the Fermi surface topology. 

The d-wave pair-pair correlation /-^(R.) is defined 
as AR.A Ri+R .), where A Ri = c RiT (c Ri+Xi + 

CRi-xi - CR i+ yx - CRi-yl)- The long range part of P<j(R) 
is a flat plateau for nonzero A, and we define P$ ve as the 
averaged value of the | R |> 2 part of P<j(R) to estimate 
the strength of SC of the system. 
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FIG. 1: Pj ve for J/t = 0.3 with several t' and t" for a hole- 
doped 12 x 12 lattice. 



FigGJshows the variational Monte Carlo (VMC) results 
of P$ ve for several t' and t" for different hole densities 5 
with J/t = 0.3 in a 12 x 12 lattice. P% ve for the t — J 
model (open circles) has the "dome-like" shape which 
is similar to the experimental results of T c versus dop- 
ing. However, this could be an artifact of the variational 
study which, we believe, overestimates the order param- 
eters and will be largely suppressed when we go beyond 
variational calculation as shown in our previous study @. 
The most suprising result is that when t' is included, P% ve 
changes dramatically. For the overdoped regime, P£ ve is 
greatly enhanced by almost one order of magnitude for 
t'/t = -0.3 and t"/t = 0.15. The SC region extends to 
5 ~ 0.4 and the peak of the superconducting dome is at 
5 ~ 0.3 and the magnitude of the maximal P$ ve is about 
2.5 times larger than for t' = t" = 0. In the undcrdoped 
region, P£ ve is almost unchanged or very slightly sup- 
pressed. Another thing to note is that beyond the value 
of t'/t = —0.3 ~ —0.4 P$ ve is no longer enhanced. 
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Fig[2l plots the maximal possible value of P°[ ve for all 
doping density as a function of t' . The maximal P% ve is 
proportional to t' in the range > t' > —0.3 ~ —0.4. 
Beyond these values pairing is no longer enhanced. Co- 
incidentally these values are about the same value of 
t'/t for mercury cuprates as estimated by Pavarini et 
aL|lO|] but much larger than what was reported in Ref. [jj. 
Among all the cuprate series, mecury cuprate maintains 
the record of having highest T c for almost a decade. 

The above result resolves the discrepancy between 
previous denisty-matrix-renormalization-group (DMRG) 
studies 0, and the band structure calculation • 
As DMRG studies were concerned with underdoped re- 
gion while the highest T c ^ max examined by the band 
structure calculation certainly depends on optimal and 
overdoped regions. 
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FIG. 2: Maximal PT" for different t' with t" = -t'/2 for 
8x8 (open circle) and 12 x 12 (full circle) lattices. 
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FIG. 3: (a)P,T e and (b) n(k = (tt,0)) vs. t' for S = 0.31 (full 
circles) and 0.083 (open circles), J/t = 0.3, t" = -t' /2. 

The different effects of t' and t" on P°[ ve between over- 
doped and underdoped regions are related to the shape 
of Fermi surface. FigOU shows the relations of P^ ve and 
n(k = (?r,0)) versus t' for S = 0.31 and 0.083. For the 
overdoped case (5 = 0.31), as —t' increases up to 0.3, 
n(k = (it, 0)) increases from less than 0.4 to larger than 
1 and P2 ve also increases sharply from less than 0.01 
to larger then 0.04. Since d— wave SC gap is largest at 
k = (7r,0), occupation of this k states by electrons en- 
hances P% ve ■ On the other hand for the underdoped 
case, n(k = (n, 0)) is almost unchanged because the oc- 
cupation n(k = (7r,0)) is already quite large (> 0.9) for 
t' = t" = and the effect of Fermi surface becomes 
unimportant. The slight suppression of P£ ve may be due 

to the destructive interference mechanism of the pair- 

I i 

hopping as suggested by Martins et aZ.|13j|. 

The decrease of P$ ve for — t' > 0.4 in the overdoped 
regime such as S = 0.31 is also likely the consequence 
of the change of the Fermi surface. n(k = (tt,0)) is al- 
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most saturated at — t' = 0.4 and remains unchanged for 
larger — t' . It is not difficult to recognize that as —t' be- 
comes much larger than t, electrons will occupy in sepa- 
rate regions around k = (±7r, 0) and k = (0, ±7r). Hence 
the Fermi surface becomes disjoint pieces. Although at 
— t'/t = 0.4, the Fermi surface is still connected but this 
tendency is already observed. The density of states be- 
gins to decrease and this is probably the reason for the 
suppression of pairing beyond —t'/t' > 0.4. 
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FIG. 4: Fermi surface of S = (a)0.19, (b)0.31 (c)0.42 (d) 0.49 
for 12 x 12 lattice, (e)optimal parameters t' v (squares) and A 
(circles). 

Fig^l shows the Fermi surface and optimal parameters 
A and t' v as a function of doping density for a 12 x 12 
lattice with J/t = 0.3, t' = -0.3 and t" = 0. In (a- 
d) the white region denotes n(k) > 1.2 and dark re- 
gion for n(k) < 0.5. The density with maximal P$ ve 
is near S opt = 0.31 as shown by the solid square in Fig^ 
The shapes of the Fermi surfaces are very different for 
S > 5 pt and S < S opt cases. Fig^Je) shows that in the 
region 0.2 < S < 0.3 although A becomes smaller, —t' v is 
still quite large and pairing is further enhanced. Doping 
beyond S > 0.3, — t' v begins to decrease quickly. This 
gives very low electron occupation at (ir, 0) as shown in 
Fig'EI c ) an d (d), then the pairing is reduced. This result 
shows that the enhancement of pairing by including t' 
is not due to larger A but from the deformation of the 
Fermi surface instead. 

Although we have emphasized the particular correla- 
tion between d-wave SC gap and electron occupation at 
k = (7r, 0) as the reason for enhancement of pairing, an- 
other familiar effect may also have played a role. It is 
well known that t' will shift the van Hove singularity [l^ 
in density of states, but it is always around k = (ir, 0). 
Hence it could be that the optimal t' v are chosen to have 
the high density of states at Fermi surface. This may be 
related to the observed extended region of flat band by 
ARPESjl^. 

The results of slave-boson MF calculation 20] in Fig[5] 
show similar behavior for the overdoped regime that in- 
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FIG. 5: Mean-field results of T c for overdoped regime with 
(/,*") = (0,0) (circles), (-0.3,0) (squares) and (-0.3,0.15) 
(diamonds). 



deed t' enhances T c . Including t" pushes the supercon- 
ducting regime to even larger doping density 6 by occupy- 
ing the momenta around (n, 0). A similar effect of t" can 
also be seen in Fig^by comparing the (t' , t") — (—0.3, 0) 
(full squares) and (—0.3,0.15) (open triangles) curves. 
Since the slave-boson method is not quite reliable quan- 
titatively in the underdoped regime, we did not show the 
values of T c in FigO However, if we do take the values 
literally, the values achieved for T c ^ max are not as greatly 
enhanced by t' as for the VMC result shown in Fig|21 
Similar results are reported by the interlayer tunneling 
model|2l|. 

Fig[f)|shows the pair-pair correlation 22J for the longest 
distance R = (1,3) for 20-site lattices obtained by the 
ED method. Pairing correlations for 2 and 4 holes are 
suppressed by t' and t" , but enhanced for the overdoped 
6- and 8-hole cases. The non-monotonic behavior of the 
overdoped cases is due to the level crossing of this system. 
If we focus on the s— like symmetry states, P% ve will vary 
monotonically in the region > t' > —0.3. The result of 
ED method is quite consistent with the variational and 
MF results that the enhancement of P$ ve by t' occurs for 
larger hole densities. 

In summary, in the optimal and overdoped regions, 
SC is greatly enhanced because of the deformation of the 
Fermi surface at these doping densities. n(k = (7r,0)) 
is enhanced by including t' and t". The occupation of 
(ir, 0) by electrons is important for the enhancement. The 
maximum enhancement of pairing correlation seems to be 



reached for -t'/t = 0.3 ~ 0.4. On the other hand, P% ve is 
not enhanced for the underdoped regime as n(k = (ir, 0)) 
is hardly affected by including t' . It is well accepted that 
the physics on the overdoped side is apparently much 
simpler in that experiments on the cuprates and theory 
for the t-J model indicate that the overdoped materials 
are very close to ordinary Fermi liquids, whence 'Fermi- 
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FIG. 6: P d (R = (1,3)) for versus t' (t" = -t' /2) for different 
2 (open circle) , 4 (full circle), and 6 (open triangle) holes in 
20 sites. The dash line shows Pd for the same symmetry of 
the 8-hole case. 



surface-based' arguments like ours are much more reliable 
than in the underdoped region. Our result shows that 
the extended t — J model naturally predicts the strong 
correlation^^ between t' jt and T c<max observed in ex- 
periments for monolayer cuprates. In addition it also in- 
dicates that further increase of t' beyond what mercury 
cuprates have most likely will not enhance T c ^ max . 

Although we have consistent results from VMC, ED 
methods and slave-boson MF calculations, the optimal 
doping density, S opt , in Fig^ is around 0.3 instead of 
0.17 obtained in the experiments and in the VMC results 
without including t' . But this is actually not a draw- 
back. As argued in References 0, 0, VMC is expected 
to overestimate the values of the variational parameters 
A which is related to the exhange energy J. Hence the 
pairing correlation is definitely much larger than that of 
a real ground state. From our previous experiences the 
5 op t seems to be always shifted to a smaller value when 
we improve the variational wave functions. Thus in the 
future work going beyond VMC calculations, we believe 
there is a better chance that we will have 5 op t closer to 
the experimental value. We also expect interlayer cou- 
pling will be important in getting the correct S op t- Now 
t' is shown to be important in enhancing pairing and it 
is also present in all high T c cuprates, the debate 0, 
about pairing robustness in the 2D t — J model without 
t' becomes somewhat irrelevant. The agreement of S opt 
between t — J VMC and experiments looks fortuitous. 

One of the consequences of our results is that the 



shape of the Fermi surface plays an important role for 
high temerature superconductors in the optimal and 
overdoped regions. Fig0] (b) and (c) show that the 
Fermi surface changes from hole-like to electron-like 
once the maximum pairing is reached and the pair- 
ing is reduced as doping increases. This is consistent 
with the ARPES results for La 2 - x Sr x CuOjM- It 
may also be related to the recent experiment ;23| which 
shows that the low-temperature Hall coefficient for the 
-Bi2>SVi.5iLao.49Cit06+<5 system exhibits a sharp change 
at the optimal doping density. Clearly this issue deserves 
more detailed study in the future. 
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